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Geometric phases in quantum mechanics play an extraordinary role in broadening our under- 
standing ol fundamental significance of geometry in nature. One of the best known examples is the 
Berry phase (M.V. Berry (1984), Proc. Royal. Soc. London A, 392:45) which naturally emerges in 
quantum adiabatic evolution. So far the applicability and measurements of the Berry phase were 
mostly limited to systems of weakly interacting quasi-particles, where interference experiments are 
feasible. Here we show how one can go beyond this limitation and observe the Berry curvature and 
hence the Berry phase in generic systems as a non-adiabatic response of physical observables to the 
rate of change of an external parameter. These results can be interpreted as a dynamical quantum 
Hall effect in a parameter space. The conventional quantum Hall effect is a particular example of 
the general relation if one views the electric field as a rate of change of the vector potential. We 
illustrate our findings by analyzing the response of interacting spin chains to a rotating magnetic 
field. We observe the quantization of this response, which term the rotational quantum Hall effect. 



I. INTRODUCTION 

In quantum mechanics the Berry phase is defined as 
a phase accumulated by the wave function during the 
adiabatic evolution around a closed path in the parame- 
ter space denoted by s. This phase can be obtained by 
integrating the Berry connection: 

^ i{^p\d^\^P) (1) 

along this path IJ. Here to shorten notations we define 
= ds^ ■ The vector s stands for an arbitrary set of 
parameters which change during the adiabatic evolution. 
These can be a coordinate of a particle, its momentum, 
magnetic or electric field, vector potential, pressure, vol- 
ume, a coupling to some external potential and so force. 
From the definition above it is clear that the Berry phase 
has a purely geometric interpretation. Together with the 
Berry connection (analogous to the vector potential 
in the parameter space) one defines the Berry curvature 
(analogous to the magnetic field): 

-F^A = d,Ax - dxA^. (2) 

From the Stokes theorem it follows that the Berry phase 
along some closed path can be found by integrating the 
Berry curvature over the area enclosed by this path. A 
non-trivial Berry phase, i.e. the phase not equal to or tt 
and a non-zero Berry curvature are in general associated 
with a broken time-reversal symmetry, because otherwise 
all wave functions can be made real [2] . 

The Berry curvature is directly related to the geometric 
tensor [3] introduced to describe manifolds of adiabati- 
cally connected wave-functions "0(5): 

XpA = (0|^5a|V) - (Vl^ |V')(V'|aA|V') (3) 

It is straightforward to check that J^^a = ^23[x^a], 
while — 3?[Xaia] defines the Riemannian metric ten- 
sor associated with the same manifold. The Riemann 



curvature associated with the above metrics as well as 
the components of the geometric tensor can serve as ob- 
servable independent measures of singularities like phase 
transitions [3]. The geometric tensor can be extended 
to mixed states by taking the statistical average of the 
geometric tensor associated with individual pure states . 

The Berry phase and related concepts found multiple 
applications in many different quantum and classical sys- 
tems In particular, it underlies the Aharonov-Bohm ef- 
fect; it shows up in transport in graphene in a quantum 
Hall regime [5]; it can emerge in photon interference of 
a circularly polarized light [3 [7]; it appears in anoma- 
lous quantum Hall effect in magnetic metals [5j, in var- 
ious magnetoelectric response functions 0, in Thouless 
pumps |10[ lllj and many other phenomena, including 
various forms of solid-state related topological phenom- 
ena like topological insulators (see Ref. [H] for a recent 
review). Most known applications of the Berry phase 
and ways to experimentally measure it rely on existence 
of free or nearly free quasi-particles which can indepen- 
dently interfere affecting transport properties. If the con- 
cept of quasi-particles is ill defined and transport exper- 
iments are not feasible, like e.g. in cold atom systems, 
the Berry phase could not be measured using traditional 
approaches. 



II. RESULTS 

The main finding of our paper is that in slowly driven 
isolated systems the Berry curvature emerges in a lin- 
ear response of physical observables to the quench ve- 
locity V = s. Specifically we show that if the quench 
velocity v is aligned along the A direction in the parame- 
ter space then measuring the response of the generalized 
force along the /i-direction: Mf^ = ~{^lJ{tf)\^f^'H\■^p{tf)) 
gives the A/i component of the Berry curvature: 

= const + J'^AWA + C(w^)- (4) 
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where we used the Einstein convention of summation over 
repeated indexes. The constant term gives the value of 
the generahzed force in the adiabatic hmit. This result 
is very remarkable since it shows that even without dissi- 
pation the leading non-adiabatic response of a quantum 
system is local in time. Indeed it is determined by the 
instantaneous velocity and the instantaneous matrix el- 
ements and the instantaneous spectrum. We sketch the 
details of the derivation of Eq. Q in the Section Meth- 
ods. Here we only point that this result is valid if either 
of the three conditions are met (i) the velocity v is turned 
on smoothly, (ii) the system is prepared initially in a state 
with a large gap, (iii) there is a weak dephasing mecha- 
nism in the system and the time of experiment is longer 
than the dephasing time. The first two conditions imply 
that the system is not excited at the initial time of the 
evolution and the last condition implies that even if there 
are small excitations in the system, they come with a ran- 
dom phase. Note that the dephasing does not have to be 
due to external noise, it can be e.g. due to averaging of 
over different experimental runs with slightly fluctuating 
durations. The equation Q applies to both gapless and 
gapped systems at either zero or at finite temperatures. 
It can be used for a single particle and for an interacting 
many-body system in a thermodynamic limit. However, 
we need to keep in mind that the Berry curvature is a 
susceptibility, in particular it can be expressed through 
the non-equal time correlation functions: ^ 



J", 



dtt{[d^'H{t),dx'Hm)o, 



(5) 



Here d^'H{t) stands for the operator d^j/K in the Heisen- 
berg representation with respect to the instantaneous 
Hamiltonian, [. . . ] stands for the commutator and the 
subindex means that the average is taken with respect 
to the adiabatically evolved state (e.g. the instantaneous 
ground state if the system is initially prepared at zero 
temperature). As usually we assume that the energies 
have infinitesimal positive imaginary parts to guarantee 
the convergence of the integral above. In low dimen- 
sional gapless systems or near phase transitions the Berry 
curvature can diverge. Then the linear response theory 
breaks down and the dependence Mfj,{v\) can become 
non-analytic. In this work we will not be concerned with 
these quite special situations. Let us point that in large 
systems if we are dealing with extensive couplings the 
Berry curvature is also extensive. Therefore possible non- 
extensive number of degeneracies of the ground state does 
not affect the result Eq. Q. However, these degenera- 
cies can affect topological protection of the integer Chern 
numbers associated with the Berry curvature which we 
discuss below. Let us also point that Eq. Q implies that 
the non-vanishing linear response coefficient indicates a 
broken time reversal symmetry in the Hamiltonian (pos- 
sibly by the coupling s^). This situation is opposite to 
that in the imaginary time dynamics, where the linear re- 
sponse for a similar protocol reveals the symmetric part 
of the geometric tensor [TSI . The latter is nonsensitive to 



the time-reversal symmetry. 

The equation Q can be interpreted as a Hall effect 
in the abstract parameter space. Indeed then v\ can be 
interpreted as a driving current in A-direction. The re- 
sponse A^^ is analogous to the electric field in the trans- 
verse direction and then is the analogue of the Hall 
resistance. Conversely by s\ we can understand the elec- 
tric field component in the x-direction (which can be 
viewed as a rate of change of the x-component of the 
vector potential) . Similarly by s^^ we can understand the 
vector potential along y-direction. Then the correspond- 
ing generalized force = —d^T-L is proportional to the 
y-component of the current. Then up to the coupling 

/h the Berry curvature is given by the Hall conduc- 
tivity and Eq. ^ reduces to the well known result (see 
Sec. VB for more details). Different variations of Eq. J4| 
also appeared earlier in other particular contexts IM| 

The proposed method of finding the Berry curvature 
has another significant advantage over traditional inter- 
ference based methods. Namely, it does not require 
strict adiabaticity, which is nearly impossible to achieve 
in large interacting systems. Contrary, the applicabil- 
ity of the linear response regime, where Eq. Q is valid, 
only requires that intensive quantities like the magneti- 
zation per unit volume remain small. Therefore there 
are no difficulties related to taking the thermodynamic 
limit with possible exceptions near singularities like crit- 
ical points where the Berry curvature can diverge. The 
relation Q gives a clear route for measuring the Berry 
curvature and consequently the Berry connection and the 
Berry phase in generic systems. A possible procedure 
can consist of evolving the system initially prepared in 
the ground state by quenching some tuning parameter s\ 
smoothly in time. Then at the time t = corresponding 
to the Hamiltonian of interest one measures the general- 
ized force corresponding to a different parameter s^: A^^^. 
Repeating these measurements at different velocities one 
can extract the slope of AA^{v\), which coincides with 
the Berry curvature. One can repeat the same sequence 
at different values of s confined within some closed area of 
the parameter space and evaluate the Berry phase using 
the Stokes theorem. The relation ^ also extends to weak 
continuous measurements. Then the Berry curvature can 
be extracted from the linear in velocity correction to A^^. 

If the parameter field s lies on an arbitrary compact 
manifold (surface) S like a sphere or a torus then the 
integral over the Berry curvature for a given state |?/;) 
forms a topological integer invariant known as the first 
Chern number: 



chi{\i^)) = — / dS^xJ^^x- 

^TT Js 



(6) 



where dS^x is an area form element in the parameter 
space. To observe the quantization It is important that 
the gap separating the ground state from the rest of exci- 
tations never closes on this surface. The regions in the pa- 
rameter space, where there are such nodes on the surface 
define the crossover between different quantized plateaus. 
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Thus Eq. Q allows one to experimentally or numerically 
map the manifold of degeneracies in the ground state 
wave function. We will illustrate this point below us- 
ing an example of an anisotropic spin chain. If in addi- 
tion the surface can be represented as an invariant closed 
manifold (such that F^x is constant on this manifold) 
then Eq. Q gives a generalization of the integer quan- 
tum Hall effect to the abstract parameter space with J-^\ 
being a product of a factor of 27r and an integer divided 
by the surface area. This setup becomes analogous to 
the Thouless pump (TUl E] with the difference that the 
quench parameter and the observable are arbitrary not 
necessarily related to the particle transport. 



III. EXAMPLES 

To illustrate how this idea works we will use several 
specific examples gradually increasing their complexity. 
First let us consider a spin one half particle in the exter- 
nal magnetic field described by the Hamiltonian: 



n 



-h ■ (7, 



(7) 



where a stands for Pauli matrices. It is well known [TB] 
that if we choose a path where hz is constant, hx = 
h± COS0, hy — h± sin0 and the angle (j) changes by 2tt 
(see Fig. [T]) then the ground state of the spin acquires 
the Berry phase: 



7 = 7r(l — COS(0)) = TT 1 



h 



(8) 





FIG. 1: Two ways of measuring the Berry phase for a spin in 
a magnetic field: (i) Traditional (left panel) where an external 
parameter adiabatically changes along a closed path, e.g. by 
varying the angle (p. Then the Berry phase can be measured 
through the interference signal between the original and the 
adiabatically evolved spins, (ii) Dynamical (right panel). In 
this setup one changes an external parameter hx linearly in 
time with some velocity Vx and measures the response of the 
magnetization my{vx)- The linear slope of the latter gives the 
Berry curvature (see Eq. Q). 



A conventional way of measuring the Berry phase 
in this setup will involve adiabatic motion of the spin 
around this path and looking into a signal sensitive to 
the interference of the rotated and the original spins. 



Using Eq. Q one can obtain the same result in a non- 
equilibrium protocol which does not directly involve any 
interference. We will illustrate this with a specific setup. 
First we prepare the system in the ground state of a 
magnetic field having a large negative x-component and 
fixed z-component at the moment t ^ t^. Then we evolve 
the magnetic field along the x-direction linearly in time: 
hx{t) = hx +Vxt. At time t = we measure the magneti- 
zation along the y-direction: rriy = {'4>(0)\ay\ijj{0)). This 
problem is simple enough so that it can be solved ana- 
lytically using the Weber functions [T^] , but the solution 
is quite involved so we solve the Schrodinger equation 
numerically instead. The result of these simulations for 
hx{to) = —99.5, hx — hx{t = 0) = 0.5 and /i^ = 1 is 
shown in Fig. [2] (solid line). The Berry curvature for 
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FIG. 2: Magnetization along the y-axis as a function of 
the quench velocity of the x-component of the magnetic field 
evaluated at — 0.5, hy — and h^ = 1. The dashed line is 
the linear response prediction given by Eq. |4] 



this system can be easily computed using e.g. explicit 
form of the ground state wave function and Eq. ([S]). In 
particular, for hy — it reads: 



J'yx — 



h^ 
2/i3 



^(hl + hl) 



3/2 



(9) 



The linear function J-yxVx (dashed line) is in perfect 
agreement with the low velocity asymptotic of the exact 
solution. From this Berry curvature we can immediately 
infer the Berry phase along the circular path by noting 
the rotational invariance of the system: 



2ir 



2(/i2 +772)3/2 



TT 1 



hz 
h 



(10) 



which is of course the correct result. If the rotational 
symmetry is broken, one would need to evaluate Tyx in 
a sufficiently dense set of parameters hx and hy enclosed 
by the closed path and evaluate the area integral over the 
Berry curvature using finite differences. 
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To get the first Chern number in this example we can 
use and hg as the external parameters keeping the to- 
tal magnetic field fixed. Then a similar quench procedure 
will result in « -T^e^e, where J^^g — l/2sin(0), which 
after integration over the spherical surface will result in 



chi — 



1 

2^ 



d0^sm{e) = 1. 



(11) 



The quantization unit here is the system's spin s — 1/2. 

Next we extend the example above to a Heisenberg 
spin chain in a magnetic field: 



N 



N-l 



2^ CTjCTj + l, 



(12) 



where N is the chain size. In this setup we will analyze 
the response to rotating magnetic field in the 6 direction 
(see Fig. [T]) fixing the magnitude h — 1: 



Kit) 



sm 



, hz{t) — cos 



, hy{t) = 0. 



(13) 

This choice of time dependence guarantees that the angu- 
lar velocity is turned on smoothly and the system is not 
excited at the initial moment of evolution. At the point of 
measurement t = tt/v the velocity of the ^-component of 
the magnetic field is exactly v. As in the previous exam- 
ple we numerically solve the time dependent Schrodinger 
equation and evaluate the magnetization along the y axis: 
TO0 = J2j i'^j) ^ ^ function of the velocity. The slope 
T^g must give the Berry curvature, which we compute as 
a function of the interaction coupling J. 

In order to check the quantization of the first Chern 
number in this system one needs to integrate the ex- 
tracted value of T^g over the sphere h = 1. However, 
in this case it is not necessary because of the rotational 
invariance of the interaction. Thus the integration will 
result in just the multiplication of J^^e by the area of the 
sphere, which is An. We thus anticipate that J-^g is quan- 
tized in units of 1/2 irrespective of the interaction cou- 
pling. This is indeed what we observe (Fig. [s]) . In the fer- 
romagnetic regime ( J > 0) the Berry curvature is equal 
to N/2, which indicates that the ground state behaves as 
a collective spin of magnitude S — N/2. In the antifer- 
romagnetic regime J is large and negative such that the 
system behaves as a spin singlet S = for even TV and as 
an effective spin 1/2 for odd A^. The transition between 
the spin singlet (for even chain) and the maximally po- 
larized state occurs through the quantized steps, which 
reflect the total value of the spin in the initial state. It is 
interesting to note that even though we numerically ex- 
tracted the slope of magnetization T^g « My (v) /v from 
a moderately small velocity v = 0.1 the accuracy of the 
quantization of plateaus is better than 0.1%. 

The example above has still one significant simplifica- 
tion coming from the fact that the magnetization com- 
mutes with the Heisenberg interaction term. Therefore 
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FIG. 3: The Berry curvature extracted from numerically 
evaluating dynamical response of the magnetization of a 
Heisenberg spin chain to the rotating magnetic field. The two 
lines represent chains of the length N = 9 and = 10. The 
quantization plateaus clearly indicate the topological charac- 
ter of the response. Numerical simulations were done by solv- 
ing time dependent Schrodinger equation with fixed velocity 
V — 0.1 (see Eqs. ( 12 1 and ( 13 1), The Berry curvature was ex- 



tracted from the transverse magnetization: 



4>0 



(v)/v. 



the time evolution of the magnetization is decoupled from 
the latter. To show that the quantization of the dynam- 
ical response holds in a more generic setup we will next 
consider a disordered (and hence nonintegrable) Heisen- 
berg chain with the Hamiltonian 



n 



N 



N-l 



(14) 



where and rjj are random variables, which for con- 
creteness are both chosen from a box distribution in the 
interval [0.75, 1.25]. We repeat the same protocol as be- 



fore, i.e. change the magnetic field according to Eq. ( 13 1, 



and extract the Berry curvature from the response of the 
magnetization to the velocity v. In Fig. [4] we show the re- 
sults of such simulations for a chain of length = 9 with 
a given realization of disorder. The slopes are extracted 
from two different velocities ?; = 0.1 and v — 0.025. 
The plot clearly shows that the quantization of the re- 
sponse persists. At higher velocity the crossovers be- 
tween the plateaus are slightly more rounded and one 
observes small fiuctuations of the numerically extracted 
T^g in the plateau regions. At smaller velocity, i.e. closer 
to the linear response regime, these fiuctuations are sup- 
pressed and we see again nearly perfect quantization. 
As a final example we analyze an anisotropic chain: 



N N-1 



3 + 1 



0.75aja|+i), (15) 



where aj^ denotes x, y components of the spins. This 
chain has only azimuthal symmetry in the xy plane. In 
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FIG. 4: Berry curvature extracted from dynamical simula- 
tions of the disordered Heisenberg chain (see Eq. (14l). The 
simulations are done for the spin chain of length A'' = 9 using 
the protocol identical to that in Fig. [3] for a particular real- 
ization of the disorder. The two lines represent two different 
velocities ?; = 0.1 and v — 0.025. Smaller velocity clearly 
improves the accuracy of quantization. 




N = 9 




FIG. 5: Berry curvature averaged over the polar angle for the 
anisotropic spin chain chain (see Eq. (15l). The simulations 



are done for spin chains of lengths N — 9 and A'' = 10 and the 
angular velocity v = 0.025. The inset shows the result for the 
even chain in the log scale. Both even and odd chain show a 
clearly quantized plateau at J > —0.25 at the value N/2. The 
even chain also shows well defined plateaus at other integer 
values. The regions between the plateaus correspond to the 
regions where the degeneracies of the ground state cross the 
surface h =const. 



additional averaging over the polar angle. Both chains 
show clear quantization of the response with F^g = N/2 
for J > —0.25. But at smaller values of J the behavior 
of the response is different. For odd chain the plateaus 
have very large fluctuations while for the even chain the 
plateaus are perfectly defined. We can understand the 
odd chain result at e.g. large negative J as coming from 
the double degeneracy of the ground state where the un- 
paired spin can be localized on the left or right edges of 
the chain (or alternatively degeneracy between symmet- 
ric and antisymmetric combinations of these spins). In a 
rotationally invariant system these states are protected 
by the symmetry and there are no transitions between 
them so the system behaves as effectively gapped. In 
an anisotropic chain this protection is lost and thus the 
ground state is not protected, hence we see no accurate 
quantization of the plateaus. For the even chain the in- 
tervals between the plateaus correspond to the regions 
where the ground state degeneracies cross the integra- 
tion surface of constant magnetic field h = 1. By de- 
forming the shape of the integration surface such that 
the jumps between the plateaus become sharp one can 
map the precise location of the degeneracies. Analysis of 
the microscopic origin of this response as well as under- 
standing which plateaus (except for the trivial ones with 
F^e/N = 0, 1/2) survive in the thermodynamic limit is 
beyond the scope of this work. 

In the examples above we focused on magnetic systems. 
But our analysis goes through if we consider e.g. dipoles 
in a time-dependent electric field or other setups. For 
isotropic systems the quantized Berry curvature reveals 
the total spin of the system. On a more fundamental level 
this quantization reveals the number of non- equivalent 
submanifolds of degeneracies of the ground state within 
the integration surface. In the analyzed examples for il- 
lustration purposes we assumed that the (^-factor, i.e. the 
coupling of spin to the magnetic field, is unity. In real ex- 
periments this quantization of the response can be used 
for e.g. a precision measurement of the g-factor like the 
ordinary quantum Hall effect is used for the precision 
measurement of /h. Quantization of the integrated 
Berry curvature can also be used for canceling effects of 
random static magnetic fields, which might affect accu- 
racy of direct measurements of the magnetization. These 
and other potential applications will be a subject of fu- 
ture investigation. 



IV. DISCUSSION AND CONCLUSIONS 



order to get the quantization of the response we thus 
need to average the Berry curvature over the polar an- 
gle of the magnetic field with respect to the z-axis. If 
we are taking a weak continuous measurement then this 
average is equivalent to the time average. The results of 
the numerical simulations for even and odd chains with 
lengths A^ = 9, 10 are plotted in Fig. [5] To obtain these 
results we used the same protocol of changing the mag- 
netic field as in Eq. ( 13 1 except that we performed an 



Our approach allows an interesting possibility of map- 
ping the Hilbert space topology through quantum dy- 
namics. By measuring the Berry curvature one can ex- 
perimentally analyze topological properties of the ground 
states, extract information about their possible degen- 
eracies and analyze phase transitions between different 
topological states. One can also use Berry curvature as a 
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probe of time reversal symmetry breaking in complicated 
systems e.g. in biology. 

In conclusion we demonstrated that the Berry curva- 
ture can be measured in generic systems, interacting or 
not, as a leading non-adiabatic response of physical ob- 
servables to quench velocity. This method does not re- 
quire stringent adiabatic conditions hard to achieve in 
large systems. While in this paper we focused on the 
quantum dynamics close to the ground state, our main 
result Eq. Q applies to the mixed states as well. We 
illustrated applicability of this method by numerically 
solving the time dependent Schrddinger equation for a 
single spin and different interacting spin chains. In all the 
cases we found quantization of the dynamical response 
in agreement with theoretical expectations. This quanti- 
zation can be interpreted as a dynamical quantum Hall 
effect. Our findings also reveal deep connections between 
quantum dynamics and equilibrium geometric properties 
associated with the adiabatically connected wave func- 
tions. We believe that these findings can be used in a 
variety of systems to experimentally analyze the valuable 
geometrical properties of the interacting systems without 
need to perform sensitive and not always feasible inter- 
ference experiments. 



V. METHODS. 

A. Sketch of the derivation of Eq. Q 

Derivation of Eq. Q is rather simple. It relies on the 
adiabatic perturbation theory, i.e. perturbation theory in 
the instantaneous basis |T71 [TH] . For linear quenches one 
finds that the transition amplitude to the n-th eigenstate 
of the final Hamiltonian is [18] 



{n\dxH\0) 



(16) 

where 8„o is the full phase difference (including the dy- 
namical and the Berry phases) between the n-th and the 
ground instantaneous eigenstates during time evolution: 



g„(y)-go(AO 

v{X') 



e„o(A) = J dx' 

(17) 

If the initial state has a large gap or if the protocol 
is designed in such a way that the initial evolution is 
adiabatic then Eq. (16) takes a particularly simple form: 



-IVx 



{n\dxH\0) 



{£n - Soy 



(18) 



The contribution of the initial term in Eq. ( 16 1 to the 



expectation value of the off-diagonal observables can be 
additionally suppressed by the fast oscillating phase Qno- 
This suppression will happen even if by magnitude this 



contribution is comparable to the final term ( [T8| . Let 
us note that in order to obtain this result from Eq. (19) 
in Ref. [T5] one needs to perform an additional phase 
transformation to undo the transformations given by Eqs. 
(6) and (14) of that work. Let us also point that Eq. ( 18 ) 
is equivalent to Eq. (2.10) in Ref. [11] for a particular 
choice of parameters. However, as we discuss above this 
result is only valid provided that the initial term in a 



more general Eq. (16) is unimportant. 

From Eq. (18) it is straightforward to derive Eq. 



- (VI - a^HlV) « (0| - d^H\0) 
{0\d^,n\n){n\dx'H\0) - fi ^ X 



{£n — £o) 



(19) 



It is easy to check that the second term in this equation 
is equivalent to Eq. 



B. Application of Eq. (Q to the integer quantum 
Hall effect. 



He we elaborate on that the quantum Hall effect can 
be understood as a particular application of Eq. ^ . This 
discussion will closely follow the first chapter of Ref. [H] . 

The Hall current in the ordinary setup can be derived 
using the adiabatic transport theory combined with the 
Kubo formula. To compute the Hall conductivity one 
usually applies the adiabatic formulas similar to those 
derived in Appendix A, to some in general interacting 
Hamiltonian T-l^Axit), Ay(t)) on a torus of size x Ly. 
The Hamiltonian H depends on the vector potential of 
the fixed external magnetic field A"^"^^ and on slowly 
varying perturbations 9 and 0, 



^,(i) = Ars + 0/L., 

Ay{t)^A'^'^^ + 9/Ly 

The current operators Ix and ly are given by 
9^' W 



(20) 
(21) 



(22) 



One imposes the external Hall voltage £ by varying the 
parameter in time. In particular from dtAy = Ey and 
—LyEy ~ £ we see that dfO = £. Thus we can view 
the Hall voltage as the rate of change of the parameter 
6. Explicitly applying the expressions for the transition 
amplitudes derived in the previous appendix, we get that 
in the linear response 



|4|m)(m|/y|0) -x^y 

{£m - £or 



(23) 



where |0) is an instantaneous ground state. The first 
term in the RHS of this equation is exactly the curvature 
of the Berry phase in the space of parameters 9 and (f>: 



((01^5^10) -(oi^a^io)) 



(24) 
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Averaging the Berry curvature over 9, cj) is equivalent 
to the evaluation of the (first) Chern character: 



chim) = 




which is an integer number. This result gives the well 
known topological quantization condition of the Hall con- 
ductivity: 

'^H = ^ch,{\Q)) (26) 
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Note added. After this work was complete we became 
aware of a recent paper [2Dj , where partially overlapping 
results were obtained. 



[1] M. Berry. Quanta! phase factors accompanying adiabatic 
changes. Proc. Royal Soc. London A, 392:45, 1984. 

[2] D. J. Griffiths. Introduction to Quantum Mechanics, Sec- 
ond Edition. Pearson, Upper Saddle River, NJ, 2005. 

[3] J. P. Provost and G. Vallee. Riemannian structure 
on manifolds of quantum states. Comm. Math. Phys., 
76:289, 1980. 

[4] L. C. Venuti and P. Zanardi. Quantum critical scaling of 
geometric tensors. Phys. Rev. Lett, 99:095701, 2007. 

[5] Y.B. Zhang, Y. W. Tan, H. L. Stormer, and P. Kim. 
Experimental observation of the quantum hall effect and 
berry's phase in graphene. Nature, 438:201, 2005. 

[6] Akira Tomita and Raymond Y. Chiao. Observation of 
berry's topological phase by use of an optical fiber. Phys. 
Rev. Lett, 57:937-940, 1986. 

[7] F. D. M. Haldane. Comment on "observation of berry's 
topological phase by use of an optical fiber" . Phys. Rev. 
Lett, 59:1788-1788, 1987. 

[8] F. D. M. Haldane. Berry curvature on the fermi sur- 
face: Anomalous hall effect as a topological fermi-liquid 
property Phys. Rev. Lett, 93:206602, 2004. 

[9] G. Ortiz, The geometry of quantum phase transitions, 
chapter in Understanding Quantum Phase Transitions, 
edited by L. Carr (CRC Press, New York, 2011). 
[10] D. J. Thouless. Quantization of particle transport. Phys. 

Rev. B, 27:6083-6087, May 1983. 
[11] Q. Niu and D. J. Thouless. Quantised adiabatic charge 
transport in the presence of substrate disorder and many- 



body interaction. J. Phys. A: Math. Gen., 17:2453, 1984. 
[12] D. Xiao , M.-C. Chang, Q. Niu, Berry phase effects on 

electronic properties, Rev. Mod. Phys. 82:1959, 2010 
[13] C. De Grandi, A. Polkovnikov, and A. Sandvik. 

Nonequilibrium quantum dynamics in imaginary time. 

Phys. Rev. A 84:224303, 2011. 
[14] M. V. Berry and J. M. Robbins, Chaotic classical and 

half classical adiabatic reactions: geometric magnetism 

and deterministic friction, Proc. R. Soc. London, 442: 

659, 1993. 

[15] J. J. Sakurai. Modern Quantum Mechanics, Revised Edi- 
tion. AddisonWesley, 2005. 

[16] N. V. Vitanov. Transition times in the landau-zener 
model. Phys. Rev. A, 59:988, 1999. 

[17] G. Rigolin, G. Ortiz, and V. H. Ponce. Beyond the quan- 
tum adiabatic approximation: Adiabatic perturbation 
theory Phys. Rev. A, 78:052508, 2008. 

[18] C. De Grandi and A. Polkovnikov. Adiabatic pertur- 
bation theory: from landau-zener problem to quenching 
through a quantum critical point. Lect. Notes in Phys., 
802:75, 2010. 

[19] Editor M. Stone. Quantum Hall Effect. World Scientific, 
1992. 

[20] J. E. Avron, M. Fraas, G. M. Graf and O. Kenneth, 
Quantum response of dephasing open systems. New J. 
Phys. 13:053042, 2011 



